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Abstract 


A  mathematical  model  for  transducer  arrays  is  developed 
which  is  not  restricted  to  fixed  velocity  distribution  trans¬ 
ducers.  The  model  is  general  enough  to  cover  most  transducers 
of  interest  for  Sonar  applications.  Its  formulation  is  also 
explicit  and  simple  enough  to  make  it  readily  usable  in  trans¬ 
ducer  and  array  design.  It  is  especially  well  suited  for 
analyzing  the  effects  of  transducer  head  flexing.  The  use  of 
the  model  is  illustrated  by  some  discussion  of  a  spherical 
array  of  longitudinal  vibrator  transducers  with  rectangular 
heads.  Numerical  calculations  are  given  here  for  one  trans¬ 
ducer  in  air,  and  for  this  case  some  results  related  to  head 
flexing  are  obtained  and  discussed. 
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Introduction 

Mathematical  models  for  transducer  array  analysis  have  usually  been 
based  on  the  assumption  of  fixed  velocity  distribution  transducers,  that 
is,  transducers  which  vibrate  with  the  same  velocity  distribution  under 
all  conditions.^  This  idealization  is  often  an  adequate  approximation 
for  underwater  sound  transducers,  but  there  are  important  cases  where 
the  velocity  distribution  of  a  transducer  depends  significantly  on  he 
medium  in  which  it  is  imnersed  or  on  circumstances  such  as  its  location 
in  an  array.  To  handle  such  cases  we  will  develop  here  a  mathematical 
model  without  the  fixed  velocity  distribution  assumption,  ’.vhich  is 
general  enough  to  cover  most  transducere  and  explicit  enough  to  be  readily 

usable  in  transducer  and  array  design.  The  model  1b  based  on  classical 

2 

modal  analysis  as  used  by  Lax  for  one  circular  plate  with  clamped  edges 
in  an  infinite,  plane,  rigid  baffle. 

We  will  develop  the  general  equations  t  the  model  in  the  first 
part  of  this  report.  They  can  be  used  to  take  head  flexing  into  con¬ 
sideration  in  transducer  and  array  design  and  to  interpret  head  flexing 

measurements  such  as  those  which  are  now  being  made  by  holographic  inter- 

3 

ferometry.  In  the  second  part  we  will  illv-trate  how  to  use  the  mod^ ' 
by  discussing  a  spherical  array  and  the  common  longitudinal  resonator 
type  of  transducer.  We  will  do  numerical  calculations  for  one  trans¬ 
ducer  in  air  and  obtain  some  results  related  to  head  flexing. 

It  must  be  remembered  that  a  general  model  is  merely  a  framework 
within  which  more  specialized  information  can  be  used  in  a  systematic  way. 


*  1  - 


Parka  Mathzmawcal  Labohatoium,  Incorporated 

ON*  KIVRR  ROAD  •  CaRLIAL*.  RAARaCHI  ARTTA 


Q1+24-SR-6 


Such  a  model  for  transducer  arrays  does  not  reduce  the  need  for  solutions 
to  all  the  specific  elastic,  acoustic  and  electromechanical  problems  which 
are  involved.  We  have  tried  to  formulate  the  general  model  in  such  a  way 
that  existing  solutions  and  information  about  specialized  problems  can 
be  easily  used.  It  is  hoped  that  this  will  make  the  model  readily  eppli- 
cable  to  practical  problems. 

The  Mathematical  Model 


General  Formulation 

Consider  an  arbitrary  array  in  which  the  radiating  surface  of  each 
transducer  is  an  elastic  body  6uch  as  a  membrane  or  plate  which  we  will 
refer  to  hereafter  as  the  head.  One  side  of  each  head  is  in  contact 
with  an  acoustic  fluid  medium,  while  the  other  side  is  isolated  from 
the  fluid  medium  but  is  i..  contact  with  the  other  parts  of  the  trans¬ 
ducer  such  as  the  electromechanical  drive  mechanism,  tie  rod  and  edge 
seals.  Thus,  acoustic  forces  act  on  one  side  of  the  head  from  all  the 
transducers  in  the  array  and  from  sources  outside  the  array.  Acting 
on  the  other  side  are  forces  from  the  electromechanical  driver  plus 
constraints  exerted  by  the  various  connections  with  the  rest  of  the 
transducer  when  the  head  moves. 

The  equation  of  motion  for  harmonic  vibrations  of  the  1th  head 
in  an  array  of  N  transducers  can  be  written 
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where 

r*  ,  *  the  position  vector  of  a  point  on  the  1th  head 

Pj  -  a  linear  differential  operator,  the  form  of 
which  depends  on  the  type  of  head 

(*)  ■  angular  frequency 

■  mas 8  per  unit  area  of  the  head 

vr  (r ;)  ■  normal  velocity  of  the  head 

f  (r*)  ■  normal  force  per  unit  area  exerted  on  the 

head  by  other  parts  of  the  transducer 

*  acoustic  pressure  on  the  head 

Some  specific  examples  of  P  are: 

p  z  -  ^  * -  y  ^  for  an  elastic  plate  (2a) 

U(\-  ff"1 ) 

Y  “  Young's  modulus 
r  •  Poisson's  ratio 
t  *  thickness 

p  *  T  V  for  a  membrane  (2b) 

T  -  Tension  per  unit  length 
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a  jH 

D-  YtK  -j^r 


for  an  elastic  bar 


(2c) 


K  ■  radius  of  gyration  of  cross  section 

(■  yVit  for  a  bar  of  rectangular  cross  section) 

These  three  examples  are  the  simplest  forms  of  the  differential  operators 

for  the  given  elastic  ouui.es.  More  general  operators  also  be  used; 

4 

for  example  ,  if  sh^ar  and  rotatory  inertia  are  included  for  an  elastic 
bar  we  have  for  harmonic  vibrations 


■  dimensionless  numerical  factor  (■  2/3  for 
rectangular  cross  section) 

yin.  ■  shear  modulus 

p  ■  density 


(2d) 


We  can  write  ihe  acoustic  pressure  in  the  form 

P(?s)  *  ^  J  J  *(r-J  G-Crj,  r.')  , 

where  each  term  of  the  sum  is  the  contribution  from  one  of  the  trans¬ 
ducers  in  the  array.  G ,  the  appropriate  Green's  function,  is  a 
convenient  representation  for  the  present,  but  the  formulation  is  not 
restricted  to  Green's  functions.  We  will  soon  express  the  acoustic 


(3) 
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terms  in  such  a  way  that  they  can  be  evaluated  by  any  method  for  calcu¬ 
lating  sound  fields  which  suits  the  specific  situation.  The  term  pr(r4-') 
represents  sound  waves  from  sources  outside  the  array  including  dif¬ 
fraction  from  the  transducers  and  baffle  of  the  array.  This  term  does 
not  depend  on  the  velocity  distribution  of  any  of  the  transducers  in 
the  array. 

Similarly,  we  can  separate  the  forces  exerted  on  each  head  by  the 
rest  of  the  transducer  into  a  part  which  depends  on  the  transducer 
velocity  distribution  and  a  part  which  does  not: 


fO'i)  JJ  „(?/)  ?/)  ds/  . 


(M 


The  first  term  is  the  blocked  force  per  unit  area  where  n(r*j  )  is  the 
electromechanical  transfer  function,  and  Ej  is  the  voltage  at  Che 
terminals  of  the  electromechanical  element  of  the  Jjth  transducer.  The 
second  term  accounts  for  all  the  mechanical  constraints  between  the 
head  and  the  other  parts  of  the  transducer,  where  ^(f^r-Jis  an  acousti¬ 
cal  transfer  impedance  function. 

When  Eqa.(3)  and  (4)  are  combined  with  Eq.(l)  we  have 

Pr  (?})•*( ?j) EJ  5^“ * ^ £(?; ,  O  J 

The  other  member  of  the  usual  pair  of  transducer  equations  expresses 
the  current  in  terms  of  Che  voltage  and  the  velocity.  For  reciprocal 
transducers  It  is 

(6) 


h  *  Ywi  Ei  *  JJ 


(5) 
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where  Y^i  is  the  blocked  electrical  admittance.  These  equations  corre¬ 
spond  to  the  general  transducer  equations  of  Foldy  and  Primakof f except 
that  they  use  velocity  and  current  as  the  independent  variables.  We 
have  also  specified  the  nature  of  the  transducer  to  a  somewhat  greater 
extent,  and  we  are  treating  arrays  rather  than  a  single  transducer. 

Note  that  the  extension  of  the  general  transducer  equations  to  arrays 
ha&  recently  been  discussed  by  Hickman,  Martin  and  Schenck.  Hickman 
has  also  given  a  general  mathematical  model  for  sonar  transducer  array 
systems  which  includes  amplifier-delay  line  networks  and  domes. 

In  the  usual  modal  analysis^  all  the  terms  in  £q.(5)  which  depend 
on  \s(r£)  would  define  the  eigenvalue  problem  to  be  solved.  In  the 
present  case  this  set  of  coupled  integro-differential  equations  is  too 
difficult  to  be  solved  directly.  Instead  we  6tart  with  the  relatively 
simple  eigenvalue  problem  formed  by  omitting  the  integrals  from  Eq.(5). 
The  solutions  of  this  problem,  which  we’ll  call  the  normal  modes  of  the 
head  and  denote  by  >  satisfy 

D'U^i)  s  ^  (r-)  (7) 

plus  boundary  conditions.  For  convenience,  we  will  let  the  aingle  sub¬ 
script  r\  represent  the  double  subscripts  which  would  be  required  for  two 
dimensional  elastic  bodies.  This  problem,  although  it  consists  only  of 
the  simplest  part  of  the  whole  problem  being  considered,  still  has  known 
solutions  for  only  a  few  elastic  bodies  with  simple  shapes  and  Simple 
boundary  conditions.  Fortunately,  many  of  these  are  useful  for  practical 
transducers.  The  in  Eq.(7)  are  the  natural  frequencies  of  the  modes 


0 


Parks  M  aimkm*tical  Lahomatomis*.  Iscr>iircK4iso 

OHS  KIVCB  HO4  0  »  Ca*U**K,  MAMirUl  HKTT* 


CA24-SR-6 


of  the  head.  These  modes  form  a  complete  6et  of  orthogonal  functions 
which  we'll  normalize  such  that 

dSj  =  M;  Snm  ,  (8) 

where  f“lj  is  the  total  mass  of  the  head.  Note  that  the  are  then 

dimensionless . 

The  normal  modes  of  the  head  can  be  used  as  the  basis  for  expanding 

the  general  solution  of  Eq.(„'.  Thus  we  write 

oo 

'<ft  =  I  '/.s  i.f-ji  .  <9) 

VUO 

where  the  mode  velocities,  v*i  ,  are  to  be  determined.  We  now  substitute 
Eq.(9)  into  Eq.(p),  multiply  ay  integrate  over  the  surface  of 

the  jth  head,  and  use  Eqs.(7)  and  (8)  to  obtain 


It  will  be  understood  that  sums  over n  or  m  include  all  the  modes, 
while  sutn3  over  L  or  j  include  all  the  transducers  in  the  array.  All 
the  integrals  on  the  right  of  Eq.(lO)  are  calculable  from  knowledge  of 
the  transducer  mechanism  and  construction,  the  array  geometry,  and  the 
operating  conditions.  This  set  of  algebraic  equations  can  be  solved 
for  the  after  truncating  to  a  finite  number  of  modes.  If  we  use  N * 

modes  we  have  N  N  equations  to  solve.  The  modes  of  each  head  are 
coupled  mechanically  and  acoustically  by  the  second  and  fourth  terms  on 
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Che  right  of  Eq.(lO).  The  fourth  term  also  couples  Che  modes  of  each 
head  to  the  modes  of  every  other  head  in  the  array. 

We  now  define  the  modal  pressure  function, 


( ?))  *  JJ\»  d  Sj^  , 

which  is  the  sound  pressure  produced  at  the  point  rj  per  unit  velocity 
amplitude  of  the  nth  mode  of  the  ich  head.  Using  Eq.(ll)  we  can  write 
the  last  integral  on  the  right  of  Eq.(lO)  as 


z.„;j 


(12) 


i  > 


O',  m 


which  is  the  mutual  radiation  impedance  between  the  nth  mode  of  the  ich 
head  and  the  mth  of  the  juh  head.  From  the  acoustic  reciprocity  theorem 
we  find 


Further,  since  is  symmetric  (antisymmetric)  if  t^n  is  symmetric 

(antisymmetric),  we  also  have 


(13) 


nv«u 

when  w  andm  refer  to  different  symmetries.  The  quantity  is 

expressed  in  terms  of  Green's  functions  in  Eq.(ll)  only  for  analytical 
convenience.  It  may  be  calculated  by  any  method  available,  and  the 
results  can  then  be  used  in  Eq.(l2)  to  calculate  the  mutual  radiation 


w 
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impedances  and  later  in  Eq.(20)  to  calculate  the  sound  pressure  in  the 
near  or  far  field. 

It  is  also  convenient  to  define 


(15) 


3  H  "<V»u  <*$  isi  * 
3  fl  fr(^)  J5J  • 


(16) 


(17) 


Then  Eqs.(lO)  and  (6)  become 

[iuM.  -l  T  X  J  „.i  *  1 1  ^  2-u]  v. 


mi  3  "mi  ,  (18) 


i  w 


h  *  X^ni  n«i  *  Ykiei  * 

Eqs.(l8)  and  (19)  are  the  main  results  of  the  general  formulation. 
From  them  the  transducer  and  array  behaviour  can  be  determined  without 
knowing  the  velocity  distributions  of  the  transducers.  It  is  necessary, 
of  course,  to  determine  the  normal  mode  functions,  the  transfer  impedance 
function,  and  the  electromechanical  transfer  function  and  to  solve  the 
associated  acoustics  problem.  To  make  these  equations  look  more  familiar 
ve  can  rewrite  them  for  the  caBe  of  one  transducer  vibrating  in  one  mode 
where  we  drop  the  t, j  subscripts  and  let  e*  m  *  o  : 
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*.  +  ho« 


(18a) 


1  »  ««ve  ♦  VbE 


(15*) 


Mw 


where 

Ctt»M  -  C  ~ZT~  +  )co  * 

Is  the  total  mechanical  impedance  of  the  transducer. 

The  sound  field  produced  by  the  array  at  any  point  ft  in  the  medium 

is 


fO0  *  X  II  *^0  Gr("R,  O 


ZZVncIIn«fFf^  -  XiL  Sit  Pnv  ^  * 


(20) 


i  n 


L  r\ 


In  this  linear  superposition  note  that  the  dependence  on  how  the  trans¬ 
ducers  are  driven  enters  only  through  the  V/n^ t  since  the  are 

independent  of  the  driving. 

The  total  time  average  acoustic  power  radiated  by  the  array  can 
be  obtained  by  integrating  the  product  of  the  pressure  and  the  norma' 
velocity  over  the  entire  surface  of  the  array: 

?  s  “TRe  15 

For  transducers  mounted  in  a  rigid  baffle  this  reduces  to  the  sum  of 
integrals  over  the  individual  radiating  surfaces: 
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?!  T  J|  W*  (  Tj)  pf  fj)  dSj 

j 


=  4-*«  I  III 

j  m  i  n 


HMki 


(21) 


The  quantity 

P.  r  4-  fa  Y  T"  X  V*  V  •  Z. 

win 

is  the  power  radiated  by  the  1th  transducer  when  all  the  others  are 
vibrating.  It  can  also  be  calculated  from  the  total  radiation  impedance 
of  each  mode. 


r*u 


•  II 


'flc 


’*N 


The  power 


*1  -  T  V»  Zmi  *  T  «♦  X  X  V,*  V  *  Z 

P  V. 


mj 


nm  j 


■mj  Mj 

is  associated  with  the  mth  mode  ot  the  1th  transducer  including  the 
coupling  of  the  mth  mode  to  all  the  other  modes  of  all  the  transducers. 
Thus  the  total  power  radiated  by  the  jtft  transducer  is 


I  P™  =  P;  • 

w 

Note  that  atK*  2*mi  are  us®ful  radiation  Impedances, 

and  that  they  are  defined  with  reference  to  the  mode  velocities. 


Summary  of  General  Procedure 

The  mathematical  model  consists  of  Eqs.(l6)  and  (19)  with  the 
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preceding  equations  from  which  quantities  such  as  Znvv,;i  can  be  calcu¬ 
lated.  The  following  steps  are  involved  in  using  these  equations; 

1)  Treat  the  radiating  surfaces  of  the  transducers  as  elastic 

bodies  and  determine  appropriate  normal  mode  functions,  >^m  )  > 

and  frequencies,  ,  from  Eq.(7)  and  the  boundary  conditions  at  the 

edges  of  the  radiating  surface. 

2)  Determine  the  acoustic  pressure  functions,  pn*t  ,  in  Eq.(ll) 

for  the  baffle  and  array  geometry  in  question  and  calculate  the  mutual 
radiation  impedances, ,  in  Eq.(l2).  For  a  receiving  problem 
calculate  the  in  Eq.(l7). 

3)  From  the  transducer  structure  determine  a  transfer  impedance 

(■A  ^  \  v 

.  „  ri  «  ri  /  >  and  an  electromechanical  transfer  function,  « C  , 
and  calculate  the  |  in  Eq.(l5)  and  the  in  Eq.(l6). 

ha)  For  transmitting  put  these  quantities  in  Eq.(l8),  limit  to  a 
finite  number  of  modes,  specify  the  transducer  voltages,  Ej  ,  and  solve 
for  the  mode  velocities,  V^j  . 

hb)  For  receiving  specify  the  electrical  termination  (for  example, 
open  circuit  with  the  1^*0  ),  use  Eq.(l9)  to  eliminate  the  from 
Eq.(l8),  and  solve  for  the  . 

5)  The  velocity  distribution  can  now  be  calculated  from  En.(9). 

6)  The  sound  pressure  in  the  near  or  far  field  can  be  found  from 
Eq. (20)  and  the  radiated  power  from  Eq.(2l), 
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7)  Eq . ( 19)  gives  the  transducer  currents  for  transmitting  or  the 
transducer  voltages  for  receiving.  It  also  gives  the  electrical  admit¬ 
tance  which  is  often  needed  in  analyzing  measurements. 

The  first  three  steps  show  how  the  analysis  has  been  divided  into 
three  parts  (elastic,  acoustic,  transducer  structure  and  mechanism)  which 
can  be  studied  separately.  The  influence  of  the  three  parts  on  each 
other  enters  in  the  fourth  step  where  the  mode  velocities  are  calculated. 
The  remaining  steps  determine  particular  aspects  of  the  transducer  and 
array  behaviour. 

In  the  second  part  of  this  report  we  will  discuss  particular  ex¬ 
amples  of  the  three  parts  of  the  problem.  We  will  also  consider  the 
simple  case  of  one  transducer  and  two  modes  and  Illustrate  solving 
Eq.(l8)  and  interpreting  the  results.  These  examples  are  intended  to 
clarify  the  steps  listed  above  and  to  help  make  the  general  model  readily 
usable  in  practical  problems. 

Examples  of  Use  of  the  Model 


Example  of  an  Array 

As  a  particular  example  consider  a  spherical  array  of  identical 
transducers  with  flat,  rectangular  radiating  surfaces  as  illustrated 
in  Fig.  1.  The  acoustic  part  of  the  analysis  of  such  an  array  has  been 
formulated  for  the  case  where  the  individual  radiating  surfaces  are  much 

0 

smaller  than  the  radius  of  t.ie  sphere  and  the  spherical  baffle  is  rigid. 
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Figure  1.  Spherical  array  of  transducer*  with  rectangular  heads 

It  is  Lhen  a  good  approximation  to  replace  the  flat  radiating  surfaces 
by  the  slightly  curved  portions  of  a  rigid  spherical  surface.  With  this 
simplification  the  modal  pressure  functions  can  be  found  by  the  classical 
method  of  expansion  in  orthogonal  functions.  The  result  is 


«o  j*. 


*  LfwCV  X  X  ^ 


/**«  W»J* 


where 


^AVrtl  *  /  J  1»f*i .fi'X*  (,i.  f;)  . 

VI 


* 
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h  is  the  spherical  Hankel  function,  k!  is  its  derivative,  Y t  is  the 
normalized  spherical  harmonic,  a.  is  the  radius  of  the  sphere,  j?w  is 
the  density  of  the  medium,  is  the  speed  of  sound  and  *  u>/ cw  . 

•4 

In  Eq.(.2)  “R^  is  referred  to  a  spherical  coordinate  system  with  origin 
at  the  center  of  the  array  and  oriented  so  that  the  point  "R  =  o.  f  ■> 

=  o  is  at  the  center  of  the  rectangular  head  of  the  i_th  transducer. 
The  normal  inodes  of  the  head,  which  would  natr rally  be  expressed  in 
rectangular  coordinates  (y^  }  ^  j)  with  origin  at  the  center  of  the  head, 
can  be  expressed  in  this  coordinate  system  by  the  transformation 


The  angles  $0  and  i f0  are  related  to  the  length,  JL  ,  and  width,  u/-, 
of  the  head  by 

&o  -  “Vlflu  , 

<fo  m  i/l  a.  . 

To  calculate  the  mutual  radiation  impedances  we  use  Eq.(22)  in 
Eq.(l2)  and  obtain 


VVU  V  > 

*■*/./  It  (23) 

V*.  v,7“  > 


It  is  necessary  to  express  ^  )  in  terms  of  the  angles  and 

ifj  of  the  coordinate  system  for  the  jth  head  which  is  rotated  with 
respect  to  the  coordinate  system  for  the  ith  head.  The  transformation 
required  can  be  written 
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S  _  f 

X  (\ 3  Y.  <*u ,  Pii ,  >  X*  .  «f  i ) 


(24) 


V*5M 


>* 

where  the  general  expression  for  “Pw.v  is  given  in  reference  9.  In 
practical  arrays  the  transducers  would  usually  be  oriented  with  the 
sides  of  their  radiating  surfaces  parallel.  In  that  case  and 


where 


<*i3 , 


,0): 


-iva4.j 


'  vv 


v)!C^i  -  y*)! 


a^v-v 


V  *-  V 


(v'  -  y  ) ) 


14  F/v'->., -V-/4  ;  v'-y^-l  ;  ,  v'*  V  , 

dy.y  3  (  — )  “vv*  »  y  <  v  , 

P  is  the  hypergeometric  function  (a  finite  polynomial),  d^-  is  the 
angular  difference  in  atimuth  and  j  is  the  angular  difference  in 
latitude  (see  Fig.  l).  Using  Eq.(cit)  in  Eq.(23)  we  find 

09  M 

I  I 


b*  ( *<«.) 

^(kl) 


l  »'v  “a .  fci .  H 


>Vmj 


(25) 


Vi  -ft 


Eq.(25)  gives  all  the  acoustical  quantities  required  for  solving 
Eq.(l8),  while  Eq.(22)  gives  the  modal  pressure  functions  needed  in 
Eq.(20).  Since  these  series  are  complicated  and  slowly  converging  for 
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fee.  >  >  I  ,  their  transformation  to  rapidly  converging  residue  series^ 

or  a  more  numerical  approach^  should  be  considered. 

There  has  been  little  done  on  the  many  specialized  radiation 

impedance  problems  which  are  involved  in  considering  non-uniform  velocity 

12 

distributions.  Mangulis  discussed  the  rocking  mode  of  a  rigid  circular 
disk  in  an  infinite,  rigid  plane  and  compared  the  results  to  the  familiar 
piston  mode.  He  defined  the  radiation  impedance  of  the  rocking  mode  in 
terms  of  moments,  but  the  result  is  consistent  with  our  general  definition 
in  Eq.(12).  In  our  notation,  with  n  s o  representing  the  piston  mode  and 
n  *  1  the  rocking  mode,  Mangulis  included  the  impedances  and 

He  did  not  discuss  -  Z  i0^,  but  Eq.(lU)  shows  that  these  vanish. 

11 

Another  example  is  Porter's  J  calculation  of  the  self  and  mutual 
radiation  impedances  of  flexural  circular  disks  in  an  infinite  rigid 
plane.  In  this  case  the  quantities  under  discussion  are  and 

with  h»o  representing  the  piston  mode  and  the  first  axisymmetric 

bending  mode.  Approximate  expressions  were  used  in  tnese  calculations 
for  the  bending  normal  mode  functions  for  supported  and  clamped  edges. 
These  boundary  conditions  prevent  the  piston  mode  from  accompanying  the 
bending  mode;  to  Chat  the  coupling  impedances  and  are  not 

required. 

Example  of  a  Transducer 

The  longitudinal  resonator  transducer  driven  by  a  piezoceramic  ring 
with  a  concentric  tie  rod  is  a  good  case  for  illustrating  the  transfer 
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impedance  function,  '  ,  and  the  electromechanical  transfer 

function,  nfrj^  .  Such  a  transducer  usually  has  waterproof  seals 
between  the  head  and  the  housing,  but  since  measurements  show  that 
these  seals  have  little  effect  on  the  vibration,  except  to  add  damping, 
we  will  neglect  them  here.  We  assume  that  the  tie  rod  contacts  the 
inside  of  the  radiating  surface  at  a  point  which  we  take  as  the  origin 
of  coordinates.  We  also  assume  that  the  ceramic  ring  is  thin  enough 
to  consider  that  it  contacts  the  radiating  surface  in  a  circle  of  radius 
r0  about  the  origin.  We  do  not  have  to  specify  the  type  of  head,  its 
shape  or  its  boundary  conditions  for  this  p  of  the  problem. 

When  the  surface  moves  the  tie  rod  exerts  a  force  on  it  at  the 
point  where  the  two  are  in  contact.  This  force  is  proportional  to  the 
velocity  of  the  surface  at  that  point.  We  can  write  the  part  of  ) 

contributed  by  the  tie  rod  as 

2.,  S(  S  cf;)  , 

where  Z.%  is  the  mechanical  Impedance  at  the  end  of  the  tie  rod  attached 
to  the  head. 

The  ceramic  ring  exerts  a  force  rt  all  points  of  the  circle  where 
it  contacts  the  head.  One  pirt  of  this  force  is  *1 ctromechenlcal ,  and 
the  corresponding  transfer  function  csn  be  written 

»  (26) 

where  Cj  is  the  radial  polar  coordinate  and  N  is  the  electromechanical 
transfar  ratio.  The  other  part  of  the  force  exerted  by  the  ceramic  at 
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any  point  of  the  circle  depends  on  the  velocity  of  the  head  at  all  other 
points  of  the  circle,  because  of  elastic  coupling  through  the  ceramic. 

A  simple  approximation  for  the  corresponding  part  of  ^  ^ r*  ^  T“j  }  la  to 
assure  that  the  force  exerted  on  the  surface  at  is  proportional  to 
the  velocity  at  Cj  .  In  this  locally  reacting  approximation  we  denote 

u  •  I 

the  contribution  to  ^  )  by 

where  Zc  is  the  mechanical  impedance  at  the  end  of  the  ceramic  attached 
to  the  head. 

The  locally  reacting  approximation  is  similar  to  considering  the 
ceramic  ring  divided  lengthwise  into  many  separate  thin  rodB  with  no 
elastic  co  ling  between  them.  It  is  not  clear  how  adequate  this  ap¬ 
proximation  is,  but.it  is  very  convenient,  because  it  avoids  treating 
another  elastic  problem  associated  with  head  flexing.  The  motion  in 
the  ceramic  is  usually  treated  as  a  one  dimensional  wave  depending  on 
position  along  the  length  of  the  ceramic.  But  when  the  head  bends  the 
motion  in  the  ceramic  also  depends  on  position  around  the  ring  except 
for  circular  heads  undergoing  axially  symmetric  bending. 

It  is  consistent  with  the  locally  reacting  approximation  to  use 

14 

expressions  of  the  form 

2i+  1  *V»L 

Z  *  p  c  A  " 

"  m  /mem  A  *L  kmL 

for  Zc  and  ,  in  this  way  other  structural  features  of  the  transducer 
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such  as  a  segmented  ceramic  stack,  a  tail  mass  and  pressure  release 

material  can  be  included  in  the  analysis.  In  this  expression  is 

the  density,  cm  is  the  longitudinal  spaed  of  sound,  A  is  the  cross 

sectional  area,  and  L  is  the  length  of  the  ceramic  or  tie  rod;  Z.^ 

is  the  mechanical  impedance  at  the  end  opposite  the  head,  end  fo  a  to/c  . 

m  r  w 

Combining  these  expressions  gives 

2,  SCf/.f^SC r.-r„).  (27) 


Using  this  result  in  Eq.(l5)  we  find 

Jnw*j  *  *  Zt 


where 


(28) 


»ir 

bnmj  JJ  'lJrjWrrro'*iSj  =  iff  f  »l  Jr*>6)  flnt(ro->^  <*6  * 

o 

The  last  step,  in  which  one  integration  is  done  in  polar  coordinates, 
holds  for  heads  of  any  shape  as  long  as  they  entirely  cover  the  end  of 
the  ceramic  ring.  Similarly,  using  Eq.(26)  in  Eq.(l6)  gives 

ITT 

vi-^r  J  '30) 

o 

Eqs (28)-(30)  give  all  the  transducer  parameters  which  are  needed  in 
Eq . ( 18) . 


Examp le  of  Normal  Mode  Functions 


Metal  plates  are  often  used  for  the  radiating  surfaces  of  underwater 
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sound  transducers.  For  circular  plates  the  three  familiar  boundary 

conditions  have  all  been  treated  to  some  extent  (clamped,*'*  supported*^ 
17 

and  free  ).  For  rectangular  plates  only  the  supported  edge  case  has 

18 

been  solved  exactly,  and  it  is  sometimes  used  as  an  approximate  model 

for  plates  fastened  in  such  a  way  that  the  edge  displacement  is  zero 

but  the  other  conditions  are  not  known.  The  rectangular  plate  with  free 

edges  is  an  important  case  for  transducer  applications,  but  the  normal 

19 

mode  functions  must  be  determined  by  approximate  methods. 

Another  boundary  condition,  different  from  the  familiar  three  and 

leading  to  a  simple,  exact  solution  for  rectangular  plates,  consists 

of  vanishing  of  the  first  and  third  derivatives  of  the  displacement  in 

the  direction  normal  to  the  edge.  This  condition  might  serve  as  an 

approximate  model  for  plates  with  edges  which  can  move,  but  are  not 

completely  free  of  forces  and  moments.  However,  it  is  quite  different 

from  the  true  free  edge  condition. 

Although  free  edge  rectangular  plates  do  not  have  truly  one  dimen- 
20 

sional  modes,  they  do  have  nearly  one  dimensional  modes  when  one  side 
is  much  longer  than  the  other.  The  known  normal  modes  for  the  bending 
vibrations  of  bars  with  free  ends  could  be  used  as  an  approximation  for 
the  latter  case.  We  will  take  these  bar  modes  as  a  specific  example 
which  is  simple  enough  to  be  discussed  here.  We  thus  Imagine  the  type 
of  transducer  discussed  in  the  previous  section  having  a  rectangular 
head  with  u r«JL  and  also  remember  that  r0  <  for  the  ceramic  ring 
to  fit  on  the  head. 
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For  a  thin  bar  in  which  shear  and  rotatory  inertia  can  be  neglected 
the  appropriate  differential  operator  is  given  by  Eq.(2c)  where  %  denotes 
position  along  the  length  of  the  bar.  Using  this  operator  in  Eq.(7)  and 
assuming  that  the  bar  is  uniform  and  has  a  rectangular  cross  section  we 


Ytl 


Tf  \n  (*)  -  ", 


Taking  the  origin  at  the  center  of  the  bar  the  boundary  conditions  for 
free  ends  are 


Hsl  -  J  In 


7-  *  ±  %  . 


The  normalization  condition  in  Eq.(8)  reduces  to 

4a 

J  nH(%)  n  (1C)  H  »  .  (3: 

-4* 

21 

Meirovitch's  discussion  of  the  normal  modes  of  the  free  bar  is 
one  of  the  most  useful,  because  it  includes  the  two  zero  frequency  modes 
representing  rigid  body  transverse  translation  and  rigid  body  rotation 
about  the  center  (note  that  Meirovitch  takes  the  origin  at  one  end  of 


the  bar) .  Any  linear  combination  of  these  two 


>e  modes  is  also 


a  normal  mode.  The  solutions  of  Eqs.(3l)  and  (32)  also  include  an  infi¬ 
nite  sequence  of  bending  modes  with  increasing  natural  frequencies. 

The  normalized  normal  mode  functions  for  the  two  rigid  body  motions 
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-j  =  0 

O 


,,ist<  de)  is  s> 


:;'.th  respect  to  the  center, 


.  ,  ..ode  cing  mode)  it-  :i symmetric . 


The  bendi,..  mot.  i  be  divideo 


. r 1 y ;  the  symmetric  ones 


i„w  =  »•, 


[«•''  K 


*  +■  vs  ~r 

cftS_n__ 


(36) 


where  is  a  root  of 
n 


t«.*K  — 2“  ■** 


1  C0S  ~f~ 

.  a  kr A  *  fe*«£ 

COS1  — S~ 


n  *  a,  S 


The  antisymmetric  bending  modes  are 


s  S  s5A  "  “ — vST 

n  st*  — 


sink— 

- s‘n  * 


w  •  3,5  (39) 


where  Va  is  a  root  of 


~  t  O.YN 


=  0 
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and 


a  tin* 


K£ 

2 


ft  •  J;  f  M  • 


(U) 


The  natural  frequencies  are  related  to  the  kn  by 

M 


w. 


(w) 


and  the  values  of  kn  are  given  by 

'.ZJL  a  t.SOb  TT 
k3jZ  «  a. so o  tt 

-  («-  ,  n  a  4,5  .  - 

Since  the  hyperbolic  functions  greatly  exceed  the  trigonometric 
functions,  the  expressions  for  a.n  simplify  and  the  exact  normal  mode 
functions  above  can  be  written  approximately  as 


\Jn)  -  \T a 


cot 


cosh 


cosK  +  COS  TC. 


r»  *  2, 4 


(36a) 


njy)  x  fT 

%Tl 


sm 


2-b-  >ifth  "X  "  k.  X 


*-x  . . "  -  - " 


n  *  5,  S’ . . . 


(39*) 


Since  the  first  normal  mode  is  a  constant  it  follows  from  the  orthog¬ 


onality  condition  '.hat 

■Va 

J  ^  (%)  d*  «  0  ,  m>0. 


(^3) 
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This  means  that  the  volume  velocity  of  all  the  rocking  and  bending  modes 
is  zero,  and  these  modes  are  relatively  poor  acoustic  radiators  when 
where  Jk  is  the  acoustic  wavelength. 


Evaluation  of  Modal  Velocities 


We  will  illustrate  solving  Eq.(ld)  for  the  modal  velocities  by 
considering  one  transducer.  Then  we  can  omit  the  i(  j  subscripts,  and 
Eq.(l8)  becomes 

V 


E  • 


m 


We  consider  the  transducer  to  consist  of  a  tie  rod  and  ceramic  ring 
attached  at  the  center  of  a  rectangular  head  which  can  be  described  by 
the  bar  normal  mode  functions.  Since  the  transducer  structure  is  sym¬ 
metric  only  symmetric'  normal  modes  will  be  excited  by  an  applied  voltage. 
Coupling  of  modes  by  the  transducer  structure  and  by  the  sound  field  will 
also  excite  only  symmetric  modes.  Using  Eq.(3^)  for  the  piston  mode  in 
Eqs.(29)  *nd  (30)  we  find 

Woe  *  l  » 

*  N  , 

Now  using  Eq.(36)  for  the  symmetric  bending  modes  in  Eq.(30)  gives 

ITT  | 

"«* vn  J 


cosW  (kmr0  eo*$)  + - <«>s  rp  cot  d) 

*#*  2. 
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where  is  the  zero  order  Bessel  function,  and  I  is  the  zero  order 
modified  Bessel  function.  When  Eq.(36)  is  used  in  Eq.(29)  the  inte¬ 
grations  can  be  carried  out,  but  the  complicated  formula  that  results 
probably  offers  no  advantage  over  direct  numerical  integration.  Note 
that  for  these  normal  mode  functions 

born  a  *W/N  , 

because  -  1  . 

Table  1  summarizes  some  of  the  numerical  results  which  are  required 
for  solving  Eq.(l8)  for  one  value  of  the  ratio  r°/jL  • 

Table  1.  Numerical  parameters  for  t°/jl  “  -212 


h  a 

Vn 

\mtO 

4  nm 

m  *H 

0  0  0 

1  1  1 

1 

-.833 

.089 

2  1.506 TT  2.28 TT* 

.184  i.21  -.833 

-.833 

.773 

- 

4  2.500TT  6.25  TT1’ 

.0082  -  .28  .089 

.089 

- 

- 

Using  the  results  in  Table  1 

in  Eq.(28)  we  have 

)«C  *  Z-A  + 

2-c  > 

m  ‘‘Hfc 

*  .173  Zc  , 

a  ^  SO  * 

-  \.z\  Z  -  .  i33  z4. 

The  fact  that  Is  much  smaller  than  r\9  or  r\t  is  understandable 
in  terror  of  the  zeros  of  the  modes  and  the  relative  size  of  the  ceramic 
ring  and  the  heed.  The  r<  »  o  mode  haa  no  zeros,  while  the  «sj  mode 
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has  zeros  at  7f  e  ±,27£i*  Since  the  ceramic  ring,  with  radius  r0  *  .21 Z  JLt 
lies  between  the  zeros  of  the  mode  the  driving  force  exerted  by  the 
ceramic  is  entirely  in  phase  with  this  mode,  and  it,  as  well  as  the  n*o 
mode,  is  strongly  excited.  The  It »  4  ,  and  higher  order,  modes  have  more 
zeros,  the  ceramic  ring  extends  over  out  of  phase  regions  of  the  head, 
and  these  modes  are  not  strongly  excited. 

Thus  it  will  be  a  reasonable  approximation  to  consider  that  only 
two  modes  exist.  Then  Eq.(44)  gives  two  equations  which  are,  using  &0  s  o> 

«.  *  (l<5) 

/  14  \ 


where  -  ^oa  and‘  Z.»o  3  21»a*  T^li8  ^  mode  Approximation  is  similar 

22 

to  the  two  part  disk  used  by  Woollett  and  Powers  to  simulate  a  flexing 
transducer  head  since  both  are  two  degree  of  freedom  approximations  to 
infinite  degree  of  freedom  systems.  To  emphasize  the  similarity  and  the 
differences  the  circuit  diagram  corresponding  to  Eqs.(45)  and  (46)  is 
shown  in  Fig.  2. 

Solving  Eqs.(45)  and  (46)  for  V0  and  gives 


y  M  u»  J  v  , 


v6*  e 


(47) 
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These  are  the  modal  velocities  that  are  produced  by  driving  the  trans¬ 
ducer  with  the  voltage  E  at  the  frequency  U)  .  All  aspects  of  the 
transducer  behaviour  can  be  derived  from  these  velocities. 


Figure  2.  Circuit  diagram  corresponding  to  Eqs.(l*5)  and  (U6) 
for  a  transducer  with  a  flexing  head. 
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This  two  mode  approximation  gives  the  first  two  normal  modes  of  the 
transducer  by  setting  E  =.  O  .  The  two  mode  frequencies  are  the  solutions 
of 

<■  J00*  -  «■  ~ZJ~  *  2**)  -  +  *  °*  (^9) 

obtained  by  equating  the  determinant  of  Eqs.(4p)  and  (46)  to  zero.  To 
solve  this  equation  for  u)  it  is  necessary  to  specify  the  impedances  Zx 
and  Zfi  and  the  radiation  impedances.  As  an  approximation  which  holds 
for  a  heavy  tail  mass  and  a  short  ceramic  stack  and  tie  rod  we  take 
and  Zc  to  be  pure  compliances: 

*2-  v  "  ~  */  w  C  t  * 

Z,  -  -  '■/“<■,  ■ 

We  also  neglect  the  radiation  impedances  which  holds  for  the  transducer 
in  air.  The  solutions  of  Eq.(49)  can  then  be  written 


Expanding  the  radical  gives 
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and 


is  the  longitudinal  resonance  of  the  transducer  (resonance  of 
the  mass  of  the  head  with  the  compliance  of  the  ceramic  and  tie  rod). 

The  first  term  of  Eq.(5l)  is  the  usual  expression  for  this  resonance 
(since  3^>co  *  ‘  ),  which  would  probably  be  used  in  a  design  neg¬ 

lecting  head  flexing.  The  other  terms  are  modifications  caused  by  head 
flexing  since  this  mode  is  coupled  to  the  bending  mode.  <*>  j  is  the 
resonance  of  the  first  symmetric  bending  mode  of  the  transducer.  Eq.(52) 
shows  that  this  resonance,  which  is  for  the  head  alone,  is  raised 
by  the  stiffening  effect  of  the  tie  rod  and  ceramic  (the  second  term), 
and  further  raised  by  coupling  with  the  other  mode  (the  additional  terms) 
The  velocity  distribution  for  these  two  bn  das  of  the  transducer  can 
be  determined  from  the  rati'' 


i.  to  hq  -*■  3..00 

}oi 


evaluated  at  to*  and  at  tw. 


(53 


Discussion  of  the  Example 

We  will  not  attempt  a  comprehensive  discussion  of  the  effects  of 
head  flexing  in  this  report.  However,  to  illustrate  the  use  of  the 
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results  in  the  previous  section  we  will  mention  some  aspects  of  the 
transducer  behaviour. 

If  we  had  analyzed  the  transducer  on  the  assumption  that  the  head 
was  rigid  we  would  have  found  the  velocity  cf  the  piston  mode  to  be 


'or 


iu)M  t 


}oe 


♦  2 


00 


and,  of  course,  there  would  be  no  other  modal  velocities  such  as  Va 
Eq.(h7)  can  be  rewritten  as 


Vn  -  V, 


or 


»  - 


I-- 


^02  + 


h„  4  ^22.  ■*“  <2-a  2, 

(}ol  +  2-02^ 


(LwM  +  K+1J(;wM'in7lI  +  J 

showing  that  the  factor  in  square  brackets  gives  the  effect  of  head 
flexing  on  the  piston  mode  velocity. 

The  effects  of  head  flexing  on  the  directivity  pattern  and  the 
radiated  power  are  of  great  practical  importance.  Under  the  rigid 
head  assumption  the  radiated  power  is 


F  *  “  V 
1  r  z  or 


■  00 


(5b) 


(55) 


(56) 


where  is  the  radiation  resistance  associated  with  .  For  the 

WO  O0 


flexing  head  Eq.(2l)  gives 


.  .a 

Vi 

v’l 

*..*T 

|v4|  Ra*  4  R«('v;va)Foa 

|V. 

\LS- 

^L|  cos  $ 

kr 

1  ^or 

R  V 

^•0  lvor 

Norl  *00 

(57) 
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Figure  3  shows  how  the  single  resonance  u9r  for  a  rigid  head  is  split 
into  a  higher  and  a  lower  resonance  Uq  and  Ui^  when  one  bending  mode 
is  considered.  It  also  shows  how  the  two  head  mode  frequencies  ioQ 
and  are  raised  to  u>t  and  when  the  head  is  attached  to  the  rest 
of  the  transducer. 

Figure  3  shows  that  for  the  flexing  head  at  low  frequency  V0  is 
very  similar  tc  Vftr,  and  Vj,  is  very  small.  At  the  frequency  ui'e  V* 
and  V2  go  through  a  maximum  together ,  because  the  modes  are  coupled, 
and  when  the  velocity  of  one  mode  is  high  it  causes  the  other  mode 
velocities  to  also  he  high.  This  is  shown  by  Eqs.(47)  and  (<46)  for 

and  V,  which  have  the  same  denominator  and  thus  the  same  maxima.  Note 

£ 

that  if  a  O  ,  corresponding  to  no  coupling  between  modes,  these 

equations  would  not  have  the  same  denominator.  If  resistance  was  included 
in  the  calculations  a  maximum  in  V0  would  still  be  accompanied  by  a 
maximum  in  ,  but  at  the  maximum  value  of  would  exceed  that 
of  Va  at  least  for  small  coupling  (jCJ>+ Z.8jL  *<  +  Ja.*.  ♦zaa,') 

as  can  be  seen  from  Eqs.(47)  and  (48).  Similarly,  both  V0  a.id  Va 
have  maxima  at  and  with  resistance  the  maximum  value  of  V2  would 
exceed  that  of  \A  . 

Q 

V0  goes  to  zero  when  the  numerator  of  Eq.(47)  vanishes,  which 
occurs  in  this  case  at  a  frequency  close  to  but  not  the  same  as 
The  numerator  of  Eq.(48)  shows  that  \ does  not  vanish  in  this  frequency 
region.  At  high  frequency  ^9/j%  approaches  n©/^  «-  J.a.  In  Fig.  3 
there  is  a  phase  change  at  each  peak  and  at  the  aero  of  V0  .  Thus  at 
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high  frequency  Ve  and  Va  are  out  of  phase. 

The  velocity  distribution  of  the  head  is 

's(y.')  *  V0in*fTO  +  ^9) 

The  velocities  at  the  center,  w*  (ei)  f  and  at  the  ends,  ,  of  the 

head  are  given  as  a  function  of  frequency  in  Fig.  4.  At  very  low 
frequency  the  motion  of  the  head  is  almost  uniform,  but  as  the  frequency 
comes  closer  to  io ^  bending  increases  and  the  velocity  at  the  ends  exceeds 
the  velocity  at  the  center  by  a  factor  of  two  or  more.  At  a  frequency 
slightly  above  Wg  the  velocity  at  the  center  goes  to  zero.  As  Fig.  3 
shows  this  is  the  frequency  where  r  l, II  and  the  two  modes  cancel 

each  other.  Above  this  frequency  the  velocities  of  the  ends  and  center 
are  out  of  phase,  and  there  is  a  node  in  the  velocity  distribution.  In 
the  vicinity  of  the  velocities  of  the  ends  and  center  have  almost 

the  same  magnitude.  At  higher  frequencies  the  velocity  at  the  center 
exceeds  that  at  the  ends  and  the  ratio  \r(o)^ w  (*/^)  approaches  -3  . 

To  show  how  the  calculations  can  be  related  to  electrical  admittance 
measurements  we  note  that  Eq.{l9)  becomes  in  this  case 

l  •  V0  +  Yb  E  .  (6°) 

If  we  take  s  io  ,  where  is  the  blocked  capacitance,  and  use 
Y  *  we  can  rewrite  Eq.(6o)  as 


CtNl  N  C<MG  N  CcNE 
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where  n*/rt  and  n*Ai  are  8iven  in  Table  1  and  the  ratios 

and  Va/ctjs|E  are  8*-vetl  by  Eqs.(h7)  and  (48).  The  magnitudes  of  these 

ratios  are  shown  in  Fig-.  2.  Note  also  that 


C.N*  '  So 

where  H  is  the  electromechanical  coupling  factor.  From  Eq.(6l)  and 
the  results  in  Fig.  2  it  is  clear  that,  when  we  neglect  all  resistance, 
we  have  large  peaks  superimposed  on  the  blocked  susceptance  curve  at 
and  If  resistance  was  included  in  the  calculations  we  would  have 

admittance  circles  in  the  vicinity  of  these  frequencies.  In  general 
we  would  have  admittance  circles  at  each  of  the  coup.ed  mode  frequencies 
of  the  transducer. 


Cone lus ion 

A  general  mathematical  model  not  restricted  by  the  fixed  velocity 
distribution  assumption  has  been  presented.  It  Bhould  be  adequate  for 
most  transducer  array  problems,  and  it  should  be  especially  convenient 
for  problems  connected  with  transducer  head  flexing.  The  application 
of  the  general  model  to  typical  situations  was  discussed  in  order  to 
further  clarify  how  it  can  be  used.  This  part  of  the  discussion  also 
focused  attention  ou  some  of  the  more  specialised  types  of  problems 
which  must  be  solved  before  general  models  can  be  used. 
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